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И ЭЛЕКТРОПРОВОДИМОСТИ ОТ ТЕМПЕРАТУРЫ 


Рассмотрено движение электропроводящей смазки в зазоре радиального подшип- 
ника при наличии электромагнитного поля с учетом зависимости вязкости и элек- 
тропроводимости от температуры. Дана оценка влияния параметров, характери- 
зующих магнитную индукцию и напряженность электрического поля, на основные 
рабочие характеристики подшипника. 
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Введение. Задача современной инженерной практики приводит к повы- 
шению требований, предъявляемых к условиям эксплуатации и эффектив- 
ности использования различных узлов современного промышленного и 
транспортного оборудования, в том числе подшипников скольжения. Од- 
ним из важных конструктивных элементов подшипников жидкостного тре- 
ния является смазочная среда. В последнее время в качестве смазочной 
среды используются жидкости, обладающие электропроводящими свойст- 
вами. Поэтому определение рациональных значений безразмерных крите- 
риев, присущих электропроводящим смазкам, обеспечивающих повышен- 
ную несущую способность при низком трении, является, на наш взгляд, 
весьма актуальным. 

Постановка задачи. Рассмотрим движение электропроводящей смазки в 
зазоре радиального подшипника при наличии электромагнитного поля с 
учетом зависимости вязкости и электропроводимости от температуры 
(рис.1). Предполагается, что подшипник неподвижен, шип вращается с по- 
стоянной угловой скоростью @. 











Рис.1. Схематическое изображение движения 
в зазоре радиального подшипника 
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Будем исходить из уравнений «тонкого слоя» для вязкой несжи- 
маемой жидкости при наличии электромагнитного поля. Эти уравнения при 
пренебрежении эффектом Холла и в случае малых значений магнитного 
числа Рейнольдса в рассматриваемом случае имеют следующий вид: 

2 
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ди’ д"? т’ 00 
/ / / 
0У, и 1 5%, 
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Здесь У’,, У›- компоненты вектора скорости; р’ - гидродинами- 


г! 


с'В’(Е'-У,В’); 


0; и’ ы пе“, д’ = ое “—®) | (1) 





ческое давление в смазочном слое; Ё = {0,0, Е" - вектор напряженности 
электрического поля; В ={В',0,0} - магнитная индукция; „- динамиче- 


ский коэффициент вязкости; о’ - электропроводимость; Т’- температура; 
’, 0 - полярные координаты; Шо - характерная вязкость; О’, - характер- 


ная электропроводимость смазки; (С - экспериментальная постоянная, вхо- 
дящая в зависимость вязкости и электропроводимости от температуры. 
Предполагается, что величины Е, В и скорость течения жидкости 
таковы, что можно пренебречь влиянием потока жидкости на приложенные 
электрическое и магнитное поля (это предположение подразумевает ма- 
лость магнитного числа Рейнольдса). При этом функции Е'(г,0) и В'(,0) 


считаем заданными, удовлетворяющие уравнениям Максвелла [1] 
аВ'=0, той’ =0. (2) 
Уравнения (2) удовлетворяются при 
_$(0) 


р 
Г 


Е'=соп5!, В’ : (3) 





здесь ф (9) - заданная функция. 


В полярной системе координат с полюсом в центре шипа уравнения 
контура шипа и подшипника приближенно можно записать в виде 


"=и, Г=и+ес0$0 , (4) 
где /, — радиус шипа; И — радиус подшипника; е — эксцентриситет. 
С учетом (2) уравнения движения (1) можно записать в виде: 

бр" _ |. ду, _1 9 >.) 











ду’ - д"? 00 м’ бр’ 
би У, 10% _ |. (5) 
о г г00 


Система уравнений (5) решается при следующих граничных усло- 
виях: 


‚=0, у,=0 при 7’ =и +есо$0 ; 


у’, =0, Ур = @, при И=и; 
2 (0) = рол). (6) 
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Перейдем к безразмерным переменным по формулам: 
у, =@ди, у, =@КУ, р =И+0’, б=И-Ю р(0)= р(2л) 





* от я * 
рее, $=4$ 
ш= ши, б'=9б, Т'’=ТХ. (7) 


Подставляя (7) в (5) и (6), будем иметь с точностью до членов 


«8 


2 
В АЕ 4 — Аф-+ №9°; а 0. (8) 
д 07’ 40 0’ 00 
ты 212 
Е 
где А = В = 50 - число Гартмана. 
п о то № 
и=0, у=0 при 1+77с050 =#(0); 
и=0, у=1 при 7 = 0; =. (9) 


Полагая, что А<1, №М<1, точное автомодельное решение задачи (8)- 
(5) будем искать в виде [2]: 


и= + 0(",0); у= = 7(-,0); "+69 =0; 





р 1%_@ & . 
= +2; и=9($); 
РЕ И 
, 
В ОЬ (10) 
о <}; © 1+7с0$0 


Подставляя (10) в (8) и (9) с точностью до членов О(пАф,), при- 


дем к следующей обыкновенной системе дифференциальных уравнений и 
граничных — к ним: 


р" = С, + Мау"; "= С, + М9 -Аф; (11) 
6 $’ =0; 6(0) =0; 6(1)=0; 
1 
У(0) =-1; 9) =0; [945 =0. (12) 


0 
Решение задачи (11)-(12) находится непосредственным интегриро- 
ванием: 


ии С; й С.е\\% вы -Ум4с | 
№ мае (2 +1) №? (е А. 
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А, _. С 
Мф 
А -С + Се“ _ Аф,е\“%® + №. 
№ (е*® —1) 
е\\% (с —49,- Се“ + Афе + Мфге\“® 


м (е*® 1) 


М + Се +С че \\е 





где С, = 


4 


1 
Так как [94 =0, то 
0 


ря ф (сме ^* — 2 Ме? + АМ, — АМ, +4 АМ — М? — Ме | (13) 


| №2 е?% — Ме — №2, +4 —2М 
Для нахождения безразмерного гидродинамического давления при- 
ходим к следующему уравнению 


12_@, С 
на № № 
Для интегрирования уравнения (14) вначале необходимо найти 4 





(14) 





как функцию от 9 $ Воспользуемся выражением скорости диссипации энер- 








ГИИ 
вв (“= т (15) 
49 Ид» ди В: 
Тогда повышенная температура будет определяться выражением 
‚ ан’ 
ат=—^, (16) 
с,О 


где О- расход; с, - теплоемкость при постоянном давлении. 


1 
О=вкб [| (с)ас. (17) 
0 


Дифференцируя обе части предпоследнего уравнения системы (1) 
и комбинируя (16)-(17) для определения и(0), приходим к следующему 


пк (18) 
40 р р 


уравнению: 


1 1 1 
где. А, = [9”ас; А [29'9'ас; А = [9?ас; К - тепловой па- 
0 0 0 
раметр, 
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ОИ, (19) 
5°с, [$ (<) 45 
0 


Интегрируя уравнение (18) от 0 до @, будем иметь 
1 
и ЕЕ ЕЕ ЕЕ: 
К(1+^,5, (0,1) + 5, (9,1) + А,5, (6,1) 
0 
40 
где 5 (0,1) = . 
Нет) о) 
После необходимых вычислений получим: 


2 п 9 
ое ла 
-П 


2 1- 1 
5, (9,1) ==> ага А в. Е : 
(1-°) 1+71_2)} 1—1 1+1с0$0 


2+1? 1-70 
5, (9,1) ато не — 
7 и 


К = 





р (20) 














ЕЕ тб й 1 | тб (21) 
2 (1-я) 1+7с0$0 2(1 —П ) (+ 7соз0) 
В качестве вязкости /{ возьмем среднее арифметическое 9 =0 и 
9 =л (в нагруженной области): 
„(+ ит) (22) 
2 
Интегрируя уравнение (14) от 0 дос учетом (22), будем иметь 
р=иС,5, (91) + иС,5, (9,1). 
Замечая, что при 9=2л арктангенсу соответствует значение Л’, 


получим необходимые для дальнейшего значения этих функций при 
9=2л: 


и 


л(2+1° 
м в 5. (2 )= 


(1 м п)” 


2л : 


Из условий (6) найдем С.: 


5, (2л,п) = 


7 2 № 
С, = 2 (6) (т -1), С, находим по формуле (13). 
(2+) 
Перейдем к аналитическим выражениям для определения основных 


рабочих характеристик подшипника: несущей способности, силы трения и 
расхода. 
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Проекции главного вектора воздействия смазок на шип определя- 
ются формулами: 


2л 2л 
Е, = к | р'зт040; Е, =-% | р’созбаб. (23) 
0 0 


Заметим, что в силу периодичности давления р’ и четности соз9 


второй интеграл соответственно равен нулю. Первый интеграл проще вы- 


числить по частям. Найдем 
2” 2п 


2л р р 
РЕЯ ] Р' (соз0) = -п,р’соз0|" +1, | ЧТ созбаб о [рсозво ы 


0 





21 ^^ 21 А 

й С с0$0 С. с0$0 

= >” | — 40+ | = И 
0 0 

Таккак | С. с0з0 = 211. || С, соз@ И ЗлпС, , 


г (1+77с056) (1-7) о (1+77с050)} (1-2) 


277С, ЗлпС, з (24) 


(и) с (-1:)” 


Сила трения и расход соответственно определяются по формулам: 
2” 2” 





то Е =при = 


х 





И = | №2 40 5 | В (25) 
0 0" с=0 б 0 
1 
О=пи' [ас . (26) 
0 


Результаты численного анализа полученных аналитических выра- 
жений для основных рабочих характеристик (24)-(26), приведенные на 
рис.2 и 3, показывают: 

1) безразмерная сила трения подшипника не имеет явно выражен- 
ного максимума и минимума. При значениях параметра, характеризующего 
магнитную индукцию № (0,0.4), имеет место резкое увеличение силы 
трения, а при № > 0.4 сила трения стабилизируется. Напряженность элек- 
трического поля А не оказывает существенного влияния на значение силы 
трения; 

2) при № (0,0.2) имеет место резкое уменьшение расхода смазки, 
при №>0.2 расход стабилизируется; электрическая напряженность А 
практически не влияет на расход. 
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Рис.2. Зависимость безразмерной си- Рис.3. Зависимость расхода ОА вы 
лы трения Раго от параметра А, ха- Ли 
Мок, раметра 4, характеризующего на- 


пряженность электрического поля, и 


рактеризующего напряженность элек- 
№, характеризующего магнитную ин- 


трического поля, и №, характеризую- 
щего магнитную индукцию дукцию 


Таким образом, полученные результаты позволяют: 

1) определить расход смазки, обеспечивающей жидкостный режим 
трения; 

2) оценить влияние безразмерных критериев Аи №, присущих элек- 
тропроводящей смазке, на основные рабочие характеристики радиального 
подшипника; 

3) найти область изменения значений магнитной индукции, где 

имеет место резкое увеличение силы трения и ее стабилизация; а также 
уменьшение расхода смазки. 
Выводы. В связи с тем, что наличие электромагнитного поля способствует 
стабилизации движения шипа, полученные результаты дают предпосылку 
решения основной задачи, связанной с устойчивостью движения шипа в 
подшипнике, работающего на электропроводящей смазке. 
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